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PACS numbers: 64.10.+h,64.70.D-,64.70.K-,71.15.Pd 


There has been a major controversy over the past seven 
years about the high-pressure melting curves of transition 
metals. Static compression (diamond-anvil cell: DAC) 
experiments up to the Mbar region give very low melting 
slopes dT m /dP , but shock-wave (SW) data reveal tran- 
sitions indicating much larger dT m /dP values. Ab initio 
calculations support the correctness of the shock data. 
In a very recent Letter [lj , Bclonoshko et al. propose a 
simple and elegant resolution of this conflict for molybde- 
num. Using ab initio calculations based on density func- 
tional theory (DFT) , they show that the high-P/high- 
T phase diagram of Mo must be more complex than 
was hitherto thought. Their calculations give convinc- 
ing evidence that there is a transition boundary between 
the normal bcc structure of Mo and a high-T phase, 
which they suggest could be fee. They propose that this 
transition was misinterpreted as melting in DAC experi- 
ments 0. In confirmation, they note that their boundary 
also explains a transition seen in the SW data Q. We 
regard Bclonoshko et a/.'s Letter as extremely important, 
but we note that it raises some puzzling questions, and 
we believe that their proposed phase diagram cannot be 
completely correct. Since the neighbors of Mo in the Ad 
series are the hep-structured Tc and Ru, we would ex- 
pect Mo to transform more readily to the hep structure 
than to the fee structure suggested by Belonoshko et al. 
Indeed, if their phase diagram were correct, it would im- 
ply that at P > 7 Mbar, there are two separate stability 
fields of fee, which would be highly unusual. Another 
puzzle is that in the DAC experiments 1 , the samples 
are reported to be in a liquid-like state on the high-T 
side of the boundary where Bclonoshko et al. propose 
that the material is crystalline. 

To check whether hep, rather than fee, could be the 
stable phase at high-P/high-T , we have calculated the 
Helmholtz and Gibbs free energies of the bcc, fee and 
hep phases of Mo, using essentially the same quasi- 
harmonic methods as used by Belonoshko et al. The 
phonon frequencies in our calculations were obtained us- 
ing the small-displacement method Q , with large su- 
percells of 216 atoms and full inclusion of thermal elec- 
tronic excitations. Our calculated bcc- fee boundary for 
P > 3.5 Mbar, where fee is vibrationally stable, agrees 
very closely with that of Bclonoshko et al. However, 
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FIG. 1: Schematic phase diagram of Mo. Dashed (green) 
curve: bcc-fcc boundary from Ref. [lj ; Solid (red) curve: 
bec-hep boundary from present calculations; melting curve is 
adapted from Ref. PJ; dotted boundaries are estimated. 

we find that hep is noticeably more stable in this re- 
gion and our bec-hep boundary is reported in Fig. 1 . 
Surprisingly, the bcc-fcc slope dT/dP is higher than the 
bec-hep one, so that there must be a bcc-fcc-hcp triple 
point at P ~ 2.9 Mbar and T ~ 4800 K. The resulting 
schematic phase diagram of Mo shown in Fig. 1 indicates 
that the transition seen in DAC experiments is indeed 
bcc-fcc. Concerning the possible liquid-like behaviour of 
the crystal above the bcc-fcc boundary, we comment on 
the harmonic vibrational instability of fee below 3.5 Mbar 
found by Belonoshko et al. Our calculations confirm this 
instability, and show that it is associated with a negative 
shear elastic constant. We find that hep is also elastically 
unstable at T = K and P < 3.5 Mbar. Above the bcc- 
fcc boundary, the crystal must therefore be highly anhar- 
monic and stabilized by entropy. The implication is that 
above this boundary the fee and hep phases have almost 
exactly the same free energy and both are on the verge 
of elastic instability. It seems highly likely that these un- 
usual features underly the controversial interpretation of 
the static compression measurements. 
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